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Abstract 

Gravitational instabilities of isothermal spheres are studied in the presence of a positive 
or negative cosmological constant, in the Newtonian limit. In gravity, the statistical 
ensembles are not equivalent. We perform the analysis both in the microcanonical and the 
canonical ensembles, for which the corresponding instabilities are known as 'gravothermal 
catastrophe' and 'isothermal collapse', respectively. In the microcanonical ensemble, no 
equilibria can be found for radii larger than a critical value, which is increasing with 
increasing cosmological constant. In contrast, in the canonical ensemble, no equilibria 
can be found for radii smaller than a critical value, which is decreasing with increasing 
cosmological constant. For a positive cosmological constant, characteristic reentrant 
behavior is observed. 

Keywords: self-gravitating gas, gravothermal instability, cosmological constant, 
reentrant phase transition 



1. Introduction 

In a seminal work 1] , Antonov described a thermodynamic instability of self-gravitating 
systems in the microcanonical ensemble, that later became known as 'gravothermal 
catastrophe' [2]. A classic review on thermodynamics and statistical mechanics of self- 
gravitating systems is the one of Padmanabhan [3] and a more recent one is written by 
Katz An extension of Antonov's instability to the canonical ensemble, named 'isother- 
mal collapse', was given by Chavanis Extended reviews on the self-gravitating gas at 
thermal equilibrium, from the statistical mechanics point of view, are given by de Vega 
& Sanchez 6, 7] and Destri & de Vega [8]. Thermodynamics of self-gravitating systems 
can be realized as the pioneering part of a, nowadays, more general, rapidly developing, 
new field of research, that is the thermodynamics of systems with long-range interactions 



In a recent letter [12| we reported on the effect of the cosmological constant to the 



Antonov's gravito-thermal instability in the microcanonical ensemble. In the present 
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work not only do we enter in the details of this analysis and report some new results, 
but also extend the analysis to the canonical ensemble |13|- The two ensembles provide 
qualitatively different results. The understanding of the cosmological constant is of great 
importance mainly due to dark energy on the one hand (positive cosmological constant) 
and AdS/CFT correspondence on the other (negative cosmological constant). For conve- 
nience we shall call a positive cosmological constant 'dS case', a negative one 'AdS case' 
and a zero one 'flat case', although we are working in the Newtonian limit, for which the 



de Sitter and anti-de Sitter spaces are more properly called 'Newton-Hooke' spaces [14 [. 

The original system P, Q under study is a spherically bounded self-gravitating gas 
in the microcanonical ensemble, i.e. the spherical boundary shell has insulating and per- 
fectly reflecting walls. The system is studied in the Newtonian limit and in the mean 
field approximation. Antonov proved that there is no global entropy maximum. Local 
entropy extrema (metastable states) exist only for ER > — O.SSSGM^ and these equilib- 
ria are stable (entropy maxima) only if po/pR < 709, where E, R are the energy of the 
system and the radius of the shell and po, pa are the density of the centre and the edge, 
respectively. Lynden-Bell and Wood Q conjectured that at the region of no equilibrium, 
that is for ER < -0.335GAP, the system would overheat and collapse. This gravother- 
mal catastrophe picture was later confirmed by numerical simulations 1^ 1^ 13, @, [3 
and has been known as 'core collapse' fi^, which plays a crucial role in the evolution 
of globular clusters. As indicated by de Vega and Sanchez [^j the collapse is a zcroth 
order phase transition, since the temperature and pressure increase discontinuously at 
the transition (the Gibbs free energy becomes discontinuous). Gravothermal catastrophe 
can also lead to the formation of supermassive black holes . 

The canonical ensemble of the system is studied by Chavanis Q . From the statistical 
mechanics point of view, the canonical ensemble in gravity cannot be properly defined 
as explained by Padmanabhan However, it can be defined formally by the use of free 



energy and can have physical realizations, as suggested in Refs. |2ll. l22j|. In Ref. |2l| . 
the interstellar medium is studied as a self-gravitating gas in thermal equilibrium with 
the microwave background. It is shown that self-gravity in the canonical ensemble can 
explain the fractal structure of interstellar medium. The very same mechanism is ap- 
plied in f22] to explain the fractal structure of the Universe, i.e. the galaxy distributions, 
assuming galaxies have reached quasi-equilibrium. Chavanis Q studied bounded isother- 
mal spheres in the canonical ensemble and found that the self-similar behavior studied in 



Refs. |21|, 122, l23[ originates in the secondary instabilities of bounded isothermal spheres 
that lead to a fragmented collapse, associated with the King's radius of the system. In 
contrast, the Jeans radius is associated with the isothermal no n- fragmented collapse that 
occurs for GM/3/R > 2.52. Our study in the canonical ensemble can be considered as a 
generalization of Chavanis' Q study in the presence of a cosmological constant. 

The growing interest on Anti-de Sitter space, due to the AdS/CFT correspondence 
and the effect of the present value of the cosmological constant on the large scale struc- 
ture of the Universe justify a stability analysis of gravitating systems in the presence of 
a cosmological constant term. In addition, in most modern cosmological models, such a 
term accompanies, the one or the other way, the evolution of the Universe from its begin- 
ning to the present. In cosmological models with a decaying vacuum energy [24, 25, 2^ 
the cosmological constant is decreasing, so that its effect could be important even for 
stellar objects in the far past ^27^. For all these reasons, we believe it is crucial to un- 
derstand the effect of (an arbitrary value of) the cosmological constant on the stability 
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of gravitational systems. In this perspective, we study the simpler possible setup, i.e. a 
spherically symmetric, Newtonian, bounded system, to gain a basic understanding of the 
effect of the cosmological constant to the stability of self-gravitating systems. 

It has been proved by de Vega and Siebert [2^ [2^ that a thermodynamic limit, 
different from the usual one, does exist for a self-gravitating gas in the presence of a cos- 
mological constant and that the mean field approximation correctly describes this limit. 
In the presence of the cosmological constant we find that, in the mean field approxima- 
tion and the Newtonian limit, the negative cosmological constant (AdS case) tends to 
destabilize the system, while the positive cosmological constant (dS case) tends to sta- 



bilize it. This result further supports recent investigations of AdS instabilities ISOl. l3lll . 
In dS case many novel features arise. The system presents a reentrant behavior 12|. A 
second critical radius, above which metastable states are restored, emerges, and in the 
canonical ensemble the system undergoes reentrant phase transition, since there appear 
two critical temperatures. Reentrant phase transitions were known to occur for statisti- 



cal systems with long-range interactions |32l l33l . |34| . |35| but not for gravitating systems 



In addition, the homogeneous solution of dS has a turning point of stability, which we 
calculate analytically, and there exist infinite non-uniform solutions in the homogeneous 
radius as well as multiple series of equilibria for any radius. 

The paper is organized as follows. In section [2] we calculate the entropy extrema, in 
section [3] we present the criteria for stability that we used in our analysis, in section |3] 
the temperature and energy of the system are calculated and in section [S] the way to 
numerically generate the series of equilibria is presented and their asymptotic behavior 
is analytically studied. In section [6] is studied the stability of the homogeneous solution 
in dS in both ensembles. The main results of the microcanonical ensemble are presented 
in section [71 where the correspondence of our 'dS case' with the Schwartzschild-dS space 
is discussed, as well. The results in the canonical ensemble are given in section |S1 

2. Entropy extremization and free energy 

Consider >> 1 identical particles (stars), bounded inside a spherical shell with 
insulating and perfectly reflecting walls. In order to calculate the entropy of the system, 
one should calculate the A^-body distribution function /Ar(ri, . . . ,r/v,pi, . . . ,pn)- This 
seems an impossible task. However, if the correlations between the particles are not 
significant and an intermediate scale where the granularity of the system can be ignored 
exists, one can work in the mean field approximation 3j using the 1-body distribution 
function f{f,p, t). This can be defined by the A^-body distribution function as 

r 

/(ri,pi,t)= fN{ri,.-.,rN,Pi,---,PN)d^r2...d^fNd^P2---d^PN (1) 

or one can think as / giving the mass dm inside a volume d^fd^p: 

dm ^ f{f,p,t)d^fd^p (2) 

We assume that all particles have mass to = 1, so that we work with the velocity v 
instead of the momentum p. Once / has been determined, the density p{r) can be found 

by 

pif) = J fd^v (3) 
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and the number of particles by 

N = I fd^rd^v (4) 



The total mass is M = Nfh. 

The question posed is which / cxtremizes the Boltzmann entropy 

S/k = - j flogfd^T (5) 

with constant energy E and constant number of particles N , where d^T = d^fd^v. For 
the case without a cosmological constant, it has been proved [l[ 01 that only spherical 
configurations maximize the entropy. For a discussion on spherical configurations in the 
presence of a cosmological constant see Ref. 13 qL We will consider only spherical, static 
distributions. In the derivation of / we follow |2|. 

Using the Lagrange's multipliers /3 = 1/kT, fi the variation condition with respect to / 
is 

SS/k - I3SE + fi6N = (6) 
The energy is E — K + U, where K is the kinetic energy 

K^\J ^'fd'r (7) 
Regarding the gravitational potential ener gy U , we have to work on the gravitational 



potential 4'{r)- In the Newtonian limit (see [37| for dynamical effects of the cosmological 



constant in the Newtonian limit), the Poisson equation in the presence of a cosmological 
constant A [s^] is 

V^(t) = 47rGp - SttGpa (8) 

where p\ = For an analytical derivation see Appendi x [Appendix B[ The validity 

of the Newtonian approximation is analytically discussed in In the Newtonian limit 
it should hold 

AR^ < 1 

while the cosmological constant is negligible if 

R 

/O > PA ^ ^R^-fT- < 1 

Rs 

where Rs — IGMjc} is the Schwartzschild radius of the system. Since for Newtonian 
systems it normally is R/Rs ^ 1, the cosmological constant is not in principle negligible. 

For spherically symmetric configurations, bounded in r G [0,i?], the potential can be 
written as 

= 07V + (9) 

with 

^^^-cT^^dV (10) 



r — r 
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PAr' (11) 



Therefore the potential energy can be written as 
Using equations ([7j) , (|T2|) the variation of energy is 



+ (13) 

Using equations (IS|) and we get 

5^/fc-/35£; + M(5iV = -y",5/(^log/+l + /?(^y + (/.^ -^^dV (14) 

So that, in order for equation ([6]) to hold for all 5f , we get 

log/ + l + /?(^y + 0^ - ^ = 0^/(r,t;)=Ae-^(^+^(''^) (15) 

for j4 = e^~^ . We get the Maxwell-Boltzmann distribution likewise 'flat case'. The 
cosmological constant enters to the equation implicitly through the potential. 
Now, the density distribution can easily be calculated 

2 

fd?v = j Ae-^e-^U'^v^Al^y e~P^ (16) 

Absorbing the constants to the initial value 0(0) of the field and the central density poi 
we finally get 

p{r) = poe-^(^M-*(o)) (^7) 
In spherical coordinates, the equation ([8]), substituting equation (jl7p . gives 

^4- i ^^^'l^ir)] = 4^Gpoe-^(^M-*(°» - SttG^a (18) 



dr \ dr 

This equation holds for r < R, where R is the radius of the bounding wall. For r > R, 
it is of course ^^-^ = and (/>, 0' should be continuous at R. Let introduce 

the dimensionless variables 

X = r^AnGpoP , y ^ P{cj^ - 0(0)) , A = ^ (19) 



Then, equation p8)) becomes 

da; \ dx 



-A[-'i-^y]=--'-^ (20) 

5' 



which we call the Emden-A equation. The initial conditions are 



y(0) = , y'{0) = (21) 

The first initial condition comes from equation (jl9p . while the second from spherical 
symmetry (gravitational field at the center is zero). We call z the value of x at R: 



z = R^/AtKJpoP (22) 

We want to generate series of equilibria, that is to find y{z) for various z , A (various 
isothermal spheres) and not just find y{x) for some A. This should be done by solving 
(EU]) for various z, A keeping N, i.e. M, constant at each case. For the 'flat' A = 
case this was very easy to perform. One could just solve Lane-Emden equation for one 
z and interpret the result y{x) as y{z) with no inconsistency, since the mass could be 
considered fixed: there is no mass scale for this system. Another difficulty that enters, 
is that varying A cannot be realized as varying since A contains po, as well, that is 
different for each equilibrium configuration. Therefore, the situation in general becomes 
rather complex. We will see in section SI how we resolved the problems by introducing 
a new parameter and constructing an appropriate computer code. 

The canonical ensemble can be studied by defining the Helmholtz free energy F = 
E~TS with use of the Boltzmann entropy ([5]). It is equivalent to work with the Massieu 
function 0, Q J = -F/T that is 

J = S-^E (23) 

It is shown by Chavanis [5}] that the maximization of J with constant T is equivalent to 
the maximization of S with constant E to first order in variations 5p. This means that 
the two ensembles have the same equilibria defined by the distribution function 

3 

/ = (^A^ ' poe-^(^-^("»e-^'3'^' (24) 

It is easy to check, performing for J the previous calculations for S", that this holds true 
in the presence of a cosmological constant, too. What is different in the two ensembles 
is the stability analysis, i.e. the second order variation of entropy and free energy as we 
will see in section [3] 



3. Criteria for stability 

We want to calculate the second order variation of entropy w.r.t. perturbations 5p. 
We follow closely Padmanabhan ^3]. Maximizing the entropy for a given distribution 
p(r) we get the Maxwell-Boltzmann distribution 

1 „2 

I{r,v) = —-—^p{r)e ^.t (25) 
(znkl ) 2 

The entropy can therefore be written as 

3Af f 
S/k = —\ogT- J plogpd'r 
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We vary p(r) and therefore (7, K through 0, T respectively keeping E and M fixed. 
Keeping E fixed gives the constraint 

5K + 5U = Q 

From this equation we derive in [Appendix C| that 

^T = -^ I d'ricjySp + ^6p6cp) + 0(3) (26) 



So that, we get 



3M / 2 \ 



6S/k + fiSM = ^i^-^j J d^r{^6p + \5p6^) - j Sr5p{l + log p) 

' ' d'^r{(l)5p + -5p54>) - / d-^r- ^' 



4r2 \ 'SM J 2 ^7 J 2p 

-p j d^rSp + 0(3) = - y" d^rSp{l + logp -p+j;) 



with the second order variation being 



(27) 



If at an equilibrium, for any perturbation Sp it is S^^^S\equii < then the entropy is 
a (local) maximum and the equilibrium is (locally) stable. If there exists one or more 
perturbations for which S^^'^SlequU > 0, then the equilibrium is unstable. Let us see how 
the sign of 6^'^'> S\equii can be deduced by an eigenvalue equation Q generated by (l?7|) . 
Since the total mass is constant, for dp should hold: 



/ d^xdp = (28) 
Jo 



Let us concentrate on spherical symmetric perturbation Sp = Sp{r) and introduce the 
mass perturbation 

q{r) = SM{r) (29) 

Then 

^P=J^^ (30) 
Anr^ dr ^ ' 

The force due to perturbed distribution ((50)' — G-^ has to be finite everywhere and 
therefore q should go like q — )■ for r — 0. This means that (/(O) = 0. Then equation 
(P5)l gives that q[K) = 0. Thus, the boundary conditions are 

(Z(O) = q{R) = (31) 
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Substituting ([29|) . (pOj) into equation ((27|) and performing several integrations by part, 
we get 



1 



SMT^ jy-^'') ^^jy''^'^y^jy^4r V8.P.2 



' ^ ^''drA\U''drJ^ + i.(-±:,^]\, (32) 



The last expression can be written as 



S^^^S^ / / drdr'q{r')K{r,r')q{r) (33) 
Jo Jo 

with 

The sign of (5'^2)5 jg therefore determined by the eigenvalues of the 'matrix' K{r,r') 

rR 

dr'K{r, r')F^{r') = ^F{r) (35) 

At this equilibrium, where there is a transition from stability to instability, it should 
be ^ = 0. If for an equilibrium is found one perturbation for which ^ > then this 
equilibrium is unstable. For an equilibrium to be (locally) stable all eigenvalues should 
be negative ^ < for all perturbations. Equation (|35|) gives 

-^i^'(''«,0.i{i5..i(3^i)}.„.-,^«(. 

with 

f-R 



I—' 


d , 


( 1 d\ 


\rr2 


'd-r^ 


\^ Anpr'^ dr J 



/ dr'^'{r')F^{r') (37) 

In ([55)) the cosmological constant enters implicitly, since: 

,, GMir) 8t:G 

4> = 2 5— PA'' 

3 

The boundary conditions of (l36|) are as given in (|3T|) 

F^{0)=F^iR)=0 (38) 

We developed an algorithm that can determine eigenvalues and eigenstates for the 
boundary value problem defined by equations (|57|) . ([55)1 . The main difficulty is that 



in V enters the unknown function F^. We resolve the problem as follows. For a given 
range of ^, the problem is solved for trial values of V, call them Vr, and then the integral 
([57]) is calculated, which gives some value V. Some value ^ is indeed an eigenvalue, only 
if V = Vt and, in this case, of course V = V = Vr- The algorithm is applied for the 
dimensionless version of p6p . namely: 

where 

and z, y, y' , B are calculated at the equilibrium. Using the algorithm we can determine 
a turning point, where ^ = 0, an instability ^ > or verify a stable branch of series of 
equilibria by checking every equilibrium point for a zero eigenvalue and for an as large 
as possible range of positive eigenvalues. We performed these tasks for every series of 
equilibria demonstrated in this paper. 

In an equilibrium for which ^ = 0, there is a transition from a stable branch to an 
unstable branch or from an unstable branch to a more unstable branch (or vice versa, 
of course). Suppose you approach the turning point from a stable series (all eigenvalues 
negative) , then at the next equilibrium point after the transition, one eigenvalue becomes 
positive. In appendix Appendix D| we show how one can determine the branch with the 



additional positive eigenvalue (instability) near a turning point, from the previous anal- 
ysis. 

Performing similar calculations for the second variation of free energy J, it is straight- 
forward to find the corresponding eigenvalue problem for the canonical ensemble: 



\Tr'^ dr \A'Kpr'^ dr 



2nF^{r)=0 (40) 



Compared to eigenvalue equation ([55]) . we see that the difference is only the absence of 
the term containing the derivative of the potential. This difference changes the onset of 
the instability for the two ensembles. 

There is a way to study stability without solving an eigenvalue problem, due to a 
classical result of Poincare [s^. In thermodynamics of self-gravitating systems, it was 
for the first time applied by Lynden-Bell & Wood [2|. Briefly, it states that a change of 
stability can only occur at a point, where two or more series of equilibria have a common 
point or where they merge into each other. As indicated by Katz [40j . practically this 
means, that in the microcanonical ensemble, the change of stability happens in the 
point where (3{E) has infinite slope, while in the canonical where (3{E) has extrema. 
Equivalently this means that E in the microcanonical or /? in the canonical ensemble, has 
an extremum with respect to some other variable (e.g. the density contrast log(po/pfl) 
or z in dimensionless variables) at the turning point. We determine the point of change 
of stability in the dimensionless variables by finding an extremum of ER/GNP in the 
microcanonical case and GMf3/R in the canonical case, with respect to the logarithm of 
the density contrast log ■ 



PR 
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4. Temperature and Energy 



At the thermodynamic equihbrium, the gas sphere has the same temperature /? ev- 
erywhere, so that it is called an isothermal sphere. We define the dimensionless inverse 
temperature ^ 

B ^ ,41, 



which can be calculated with use of the dimensionless variables (jl9|) . by integrating the 
Poisson with A equation (|18p : 



dr 
R 



R^ 

= GAI -SttGpa — 
; B{z) ^ zy' + ]-Xz^ 



The kinetic energy per particle is 



GM 



z dy 
^ di 



inG 



-PA 



AnGpofi 



(42) 



K Jf^d^v /, 



■^d^v 



N 



Since M — Nrh = N we get 



-dH 



3 

2p 



K = 



3M 



(43) 



To calculate the total energy, we will use the Virial theorem, so that to avoid performing 
one more numerical integration to calculate the Newtonian potential energy and therefore 
improving the computer's performance. However, in the Appendix [Appendix E we 



perform a straightforward calculation of the energy in order to numerically cross-check 
the two expressions for some cases, who are proven to give identical results. This confirms 
the fact that bounded isothermal spheres are virialized in the presence of the cosmological 
constant. 

The Virial theorem for a discrete collection of matter: 



2 < if >= 



N 

fc=i 



rk > 



can be generalized in our case as 



2K 



j p{y(j)-r)d^r + -iPV 



(44) 



where V — ^t^R^ is the volume of the shell and P the pressure exerted on the gas by the 
walls. The term 3PV arises simply since F ■ R — P ■ AnR? ■ R = 'iPV . The contribution 
of the Newtonian potential to the right-hand side is just the Newtonian potential energy 
?7jv, while for the cosmological potential (/)a we get the term 



p(y(j)A ■ r)d^f 



SttG 
P—^PAr ■ r 



d^r^ 
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Substituting everything in equation (I44p we get the foUowing form of the Virial theorem: 



2K + Un- 2Ua = 3PV (45) 

Ehminating J/at we get 

E = 3PV -K + Wa (46) 
Using (jl9p we write the potential Ua in the dimensionless variables 



'^-PAf 47rr dr = / poe ^47r— — - — —x —— 

3 3 po 2 J (ATrGpofiy z/R 

^ e^^x'^dx — — — ——-7r — - — f e~^x^dx 



6z Gl3^ /R J 6z {GMP/RY R 

^ Ua = / e-^x-^dx (47) 



GAP 6B^z _ 

Let us calculate the term SPV in dimensionless variables 

p_ P _ Poe-y _ z^ e-y _ z^e-y 1 GAP R _ z^e-y 

We define the dimensionless energy 

We can calculate Q from the virial equation (|46l) 

^^Qi^)^—^- — -—^Je yxdx (49) 

For simplicity let us focus in the microcanonical ensemble for the moment. The self- 
gravitating gas is characterized by two instabilities [3]; a 'strong' instability that is associ- 
ated with complete absent of equilibria and a 'weak' instability which refers to equilibria 
that are unstable, i.e. entropy, although it is an extremum, is not a local maximum. 
The equilibrium point at which the weak instability sets in, called the turning point, is 
therefore the point at which the second variation of entropy becomes zero, as calculated 
in the previous section. However, Poincare's theorem insures [s^ \^ that this point is 
the same with the marginal point of the strong instability. In Figure [U where the series 
of equilibria f3 = I3{E) is drawn for the flat A = case, this is point B. The strong 
instability corresponds to the region at the left of the vertical dashed line that crosses 
B. The weak instability refers to the branch where S is the focal point of the spi- 
ral. Thus, point B is simultaneously the turning point of stability (from stable branch 
AB to the unstable branch BS) and the marginal point of the strong instability. This 
str ong instability leads to a core-halo structure as verified by monte-carlo simulations 
d, m, [13, [3| and is associated with a collapse phase transition 0, H , while 
the weak instability leads to a fractal structure and is associated with a fragmented col- 
lapse, called a clumping phase transition 0, 0, H, HI]. This fractal structure is due 
to the secondary instabilities that set in at points 5i, 5*2, etc. All of the above hold in 
the canonical ensemble, as well, where the axes in Figure [T] should be interchanged and 
entropy should be replaced with free energy. 
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1.6- 

1.4- I 

1.2 1 ' 1 ' ' ' ' ' 

-0.5 -0.4 -0.3 -0.2 -0.1 0.1 

ER/GM^ 

Figure 1: The series of equilibria /3 = I3{E) for A = 0. In the microcanonical ensemble, AB is the stable 
branch and B is the turning point of stability. As we can see, it is also the marginal point for which 
equilibrium states do exist. 



5. The series of equilibria and asymptotic behavior 

We want to solve the Emden-A equation ([20)) for various pA keeping M constant and 
for various isothermal spheres with radius R. The cosmological constant introduces a 
mass scale 

4 c. 
Ma = -ttR^ ■ PA 



We define the dimensionless mass 



= (50) 



2Ma 2pA 

where p — M/ {^ttR^) is the mean density of matter. This gives 



Stt R^pa 

that is 



m = -^-^ — (51) 



po 1 3GM/3 3B 



2pa4ttGpo(3R'^ R "" Az2 

Equation (j52p implies that in order to keep m fixed, A has to be different at each z. 
We developed an algorithm that solves the Emden-A equation keeping the quantity m 
constant; for each z, some A values are iterated until the requested value of to, calculated 
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Figure 2: For pA > (dS) there exist various series of equilibria for some fixed pA- 



by equation (l52l) . is found for a relative tolerance predetermined by the user. For the 
various plots in this paper, we used relative tolerance between 10~^— 10""'^^ depending 
on the needs of each case. From equation ([^T|) . it is evident that solving for various fixed 
m can be interpreted as varying p\ and/or R for a fixed M . Therefore we can determine 
how various quantities change with respect to p\ by solving for various m. 

For the AdS case (pA < 0) we find that for each pA only one series of equilibria exists, 
likewise fiat (pA = 0) case. However, for dS (pA > 0), we find that for a fixed pA, there 
exist more than one series of equilibria. This is evident in Figure [5] where m is plotted 
w.r.t. A, z. We see that the intersection of a plane m — const with the m-surface defines 
various different curves in the (A, z) space. 

The Emden equation (j20p for p\ — 0, i.e. A = 0, is well known to have an exact but 
singular solution 

Vs = log y (53) 

with infinite density at the origin since e^^" — It can be shown [2] that for z — > oo 
the series of equilibria approach the singular solution. Therefore, the singular solution 
in the fiat case is the focus point of the central spiral /3{E) in Figure [31 We can see 
in this picture that in dS and AdS cases there seems to exist an equivalent to the fiat 
singular solution, which can be identified with the focal points of the spirals (note that 
in dS, not all series of equilibria do form a spiral /3(E)). Unfortunately, no analytical 
solution of equation (j20p is known. However, following the fiat case paradigm we can 
determine asymptotically the equivalent of the singular solution in dS and AdS. Applying 
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ER/GM^ 
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Figure 3: The temperature versus the energy for some negative ('AdS'), zero ('Flat') and positive ('dS') 
cosmological constant. The mass At and the radius R are held constant. 

the transformations 01 

C = log2; , u = -y + 2C (54) 
the Emden-A equation becomes 

^ + -+e«-2-Ae^C = o (55) 

This is the equation of an oscillator in the potential V{u) = e" — 2w with external force 
Fext = Ae'^^ — u' where prime denotes differentiation w.r.t. C,. The external force is 
damping in AdS case (A < 0) and forced, damping in dS (A > 0). The potential has 
a minimum uq = log 2. In dS case, if the external force is damping, i.e. the term —u' 
dominates, then for ^ — )■ cx), we have u — >■ uq. Therefore, in this case, we can make an 
expansion of u about uq. The damping dominates if: 

Ae^C « u' ^ Xz^ « zy' ^ Xz^ « (56) 

where y' denotes differentiation w.r.t. z and we have used equation (I52p . This limit 
means that we are considering very small pA, just about the flat case, since Xz^ — > has 
to be taken together with z — > cxi. The two limits are consistent with each other, since 
by equation (|52[) we see that as z — >■ oo it must be A ^ so as Az^ to remain finite in 
order for m to be finite (and B is bounded for the equilibria we are considering). 
Making the transformation 

u = log2 + ui 
14 



equation ((55l) becomes 

Provided condition (j56p holds and for C — )■ cx), ui is small and we can expand the 
exponential keeping the first two terms, to get: 

d?'U^ du^ OA 

The solution of the corresponding homogeneous equation is well known to be 

Uh^Ae 2cos( — C + (5) 

and one solution of the full non-homogeneous equation can be found easily to be Up — 
je^^- We get in total 

ui = Ae~i cos(^C + 6) + ^e^^ 
2 8 

Since 

u = ui+ log 2 => — y + 2C = ui + log 2 ^ y = log ui 

the asymptotic behavior of the Emden-A equation for z — > oo and for small is 

ya = log— rC0s(— logz + d) - -z (57) 

The density p is given by the exponential e~^. We get the asymptotic behavior 

g-Va ^ ^^^z^ ^Az-^'^ cosi^^ log Z+S) 

:^ ^e^^' 1^1 + 4 cos(^ log z + ^)^ (58) 

We see that the equivalence to the singular solution (15^ of the fiat case, is given in dS 
and AdS by the 'asymptotic singular' solution e^-^^ = -^e^^ that is 

y^s = logy-^z^ (59) 

This solution uas, is easy to check that indeed satisfies the Emden-A equation (j20D to 
first order in Az^. 

By equations and we can calculate the temperature and energy of the asymptotic 
singular solution to be: 

GM(3as ^ 4m REas ^ 2e^^l^^ J 1 9_ 

R 4m - 1 ' GAf2 52 25 mS 40m2 ^ ^ 

Using this equation for Bas-, the condition gives the values of m for which the 
asymptotic singular solution is valid 

m < -6.98 and m > 5.23 (61) 
15 
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Figure 4: Graphical solution of the problem ta,n{z) = z. 

Recall that m — p/2p\ is used to vary the cosmological constant. So that m — )■ cx) 
corresponds to the flat case, m < to AdS and m > to dS. We find that the singular 
point determined analytically by equation (j60p does indeed coincide with the numerically 
determined focal point of the spiral f3{E) for various arbitrary values m in the allowed 
range (|6ip in dS and AdS. Even for values of m outside the allowed range (|6ip (like 
TO = 2 and to = —2 of Figure [3]), the analytical expression gives a result very close to 
the numerical calculation. 

6. Homogeneous solution in dS 

For a positive cosmological constant, the Emden-A equation has solution with a uni- 
form density (for r < R), we call homogeneous solution. For p — const, the Poisson 
equation gives (pN — ^^pr^, so that the full potential with the cosmological constant is 




The Poisson with A equation ([5]) gives cj) — 

P = 



2pAy + 0(0) (62) 
0(0) = const, for 

2pA 



which of course is consistent with equation (j62p . This is the homogeneous solution with 
dimensionless temperature 

Bh = (63) 
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Figure 5: The energy versus the density contrast at the homogeneous radius R = Ru for fixed M and 
positive PA- Unlike flat case, this plot cannot be realized as keeping E constant and varying R, because m 
is held fixed. We see two distinct series of equilibria that are stable for po < PR ^^nd unstable for po > PR- 
The change of stability happens at po = PR, which corresponds to specific homogeneous solutions, that 
indeed have a zero eigenvalue at these specific points {ER/GM^ = -0.677 and ER/GM'^ = -0.825). 



and dimensionless energy 

Since (/)' — the homogeneous solution has the same turning points in the two ensembles, 
because the eigenvalue equations (j36|) and (|40l) are the same in this case. Therefore, 
although the following analysis is done in terms of the microcanonical ensemble, the 
results hold for the canonical ensemble, as well. So that all turning points for the various 
solutions presented in this section are identical in the two ensembles. 

The radius Rh of the homogeneous solution is independent of energy and temperature 
and is given for a fixed mass M and cosmological constant pA by equation: 



The homogeneous solution resembles the Einstein's static universe (|4l|,|42| has found the 
Einstein's static universe to be a local entropy maximum among other possible universes) 
in the Newtonian limit. However, to be more precise there are many homogeneous con- 
figurations (uniform density) with different temperatures, from which only one is literally 
static, that is the solution with T = (/3 — > oo). The question is whether all, some or 
none of these solutions are stable. The static solution is surely thermodynamically un- 
stable since it corresponds to only one microstate, therefore it has the minimum possible 
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entropy, i.e. zero entropy. On the other hand the solution /3 = 0, which is allowed from 
(|63)) . is stable since it behaves as an ordinary gas. More formally, equation ([32]) can be 
written as 




/ drc^'q] +^ dr\- dr-^ (66) 



which gives 

< for ah perturbations q, for /3 — (T — > oo). In dimensionless 
variables equation (|66| reads 



MS'^^^S = -l(^£ dxy'(^ + ^Bz £ dx^{q^ - q'^e^) (67) 

which gives S^^^^S < for S = 0, as well. Therefore, there exists a point of change 
of stability somewhere between B = and B ^ oo. From equation (IMl) we see that 
the energy has not an extremum. However, Poincare's criterion does not exclude the 
possibility of having a change of stability at a point other than an energy extremum f ioj . 
This is our case. The differential equation for the homogeneous solution y = y' — 
(and m = X ~ 1) and for ^ = becomes 

F" --F' + F = 

X 

which for F{0) ~ has solution 

F{x) — c{—xcos{x) + sin(x)) 

We want the smallest z that satisfies the second boundary condition F{z) = 0, that is 
the change of stability happens at this z that is a solution of the equation 

tan(z) = z 

The solution, call it zq, can be found graphically (see Figure H]) to obtain 

zo ~ 4.4934 

which corresponds to Bq ~ 6.73 from equation Therefore, the homogeneous solution 
is stable for temperature T > Tq and unstable for T < Tq, with 

^° " 6J3i?^ 

As we can see in Figure [5] there are infinite series of solutions for the homoge- 
neous radius Rh that corresponds to to = 1 (the homogeneous solution corresponds 
to m = A = 1 ) . In Figure [5] is drawn the dimensionless energy Q — ^"pi versus the 
density contrast \og{po/ pji) for two of these series. We see that there are solutions with 
Po < Pb. that continuously turn to solutions with po > Pr- At the point po = PR, which 
corresponds to a homogeneous solution, there occurs a change of stability for the two 
solutions. That is because the corresponding energies Qq = —0.6771, Qi = —0.8246 
correspond to the two first zero eigenvalues zq = 4.4934, zi = 7.7251 of the homogeneous 
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Figure 6: In dS, the dimensionless energy versus the density contrast for fixed M, R, pA in the micro- 
canonical ensemble. Unlike fiat case, this plot cannot be interpreted with varying R and fixed E, because 
at each point, m is held constant. Only the upper series (1 and 4) correspond to p(r) monotonically 
changing. Distinct series corresponding to more negative energies have more extrema of p{r). At points 
B, except B4, an instability sets in. Curves AgBs and A5B5 are already unstable solutions. 



solution as can easily be verified by equation (I64p . We numerically determined the un- 
stable branch to be the one with po > PR- One is forced to conjecture that this pattern 
of different series, for R = Rh, is infinitely continued (as indicated by Figure [2|) for lower 
and lower energies. The change of stability at po — pa should always correspond to 
Q > —0.9 as indicated by equation 



7. Microcanonical ensemble 

In this section, we review our results of Ref. Let investigate the solutions for 

R < Rh and R > Rh in the microcanonical ensemble. The energy versus the density 
contrast is plotted in Figured The upper series labelled 1 (i? < Rh) and 4 (i? > Rh) 
correspond to monotonically changing density. For series 1 it is decreasing (po > Pr), 
while for 4 increasing {po < pn). Series 1 suffers a change of stability at point Bi (stable 
branch is AiBi), while series 4 is stable everywhere and does not suffer any change of 
stability. This is proved as follows: for this series the limit E' — > 00 does exist, which 
corresponds to /? = 0. By equation (|66)) for /3 = we get < 0. In addition, the 

energy does not have an extremum (where a transition to instability could occur) and 
(to be sure) the whole series is numerically checked at each point for a zero eigenvalue. 
No one is found. Every such solution (series 4) corresponds to configurations somewhat 
hollow at the center with matter concentrated mainly at the edge. The next series at 
more negative energies have one density extremum and at the next, one more is added 
and so on. At points B, except ^4 an instability sets in. For series A2B2 and AqBq 
we have strong numerical evidence that are stable, while series 3 and 5 are found to be 
unstable. Series like A2B2 and AqBq correspond to solutions diluted at the center with 
periodic condensations away from the center. One would normally expect this pattern 
in Figure [6] to continue as one finds series with more and more negative energy. 
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7.1. Critical quantities 

A general result is that as the cosmological constant increases, isothermal spheres 
exist at even lower temperatures and energies. This can be seen in Figure [31 where the 
classic spiral /?(£') is drawn for some positive and some negative cosmological constant. 
One could say that AdS destabilizes, while dS stabilizes the system. 

In the flat case, isothermal spheres exist only for 

E R> -0.335GAf2 (69) 

This means that for a fixed radius, there exists a minimum critical energy Ecr — 
~0.335GM^ / R down to which, equilibria do exist. We want to determine how this 
critical energy changes with respect to pA- The answer lies in Figure [T] The critical en- 
ergy decreases for increasing cosmological constant. In AdS the critical energy becomes 
positive for p\ < — 4.2p, where p is the mean density. 

Assuming the energy E and radius R are fixed at values that respect in fiat case, 
does not guarantee that the equilibrium is stable. This depends on the density contrast 
Pq/pr, that is the ratio of the central density versus the edge density. The situation is 
similar in the presence of pA- For p\ = it is {po/pR)cr = 709 p-j with the unstable 
branch being the one with po/pR > 709. In Figure [8] we see how this number changes 
w.r.t. the cosmological constant. The critical density contrast decreases for increasing 
cosmological constant. The instability occurs in AdS at more condensed configurations, 
while in dS at less condensed configurations. 

Assume that the energy is negative and fixed at some value E in the fiat case. Then, 
inequality (|69| shows that there is a maximum critical radius Ra = (— 0.335/i?)GAf^ 
that constrains the existence of an equilibrium. For R > Ra there are no equilibria. In 
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Figure 8: The critical density contrast versus pA for fixed M , R in the microcanonical ensemble, where p 
is the mean density of matter. An instability sets in for po / pji > (po / pii)cr (except when (po / pR)cr < 1) 
where a clumping phase transition occurs. 



AdS this radius decreases as pA attains more negative values. In dS this radius increases 
as the cosmological constant increases and in addition there appears a second critical ra- 
dius, we call RIA^ that constrains the existence of equilibria from below. At Ra and Ria 
a collapse phase transition takes place. That is, no equilibria exist for Ra < R < Ria 
as can be seen in Figure IH] and the system lies in a collapsed phase. This is a typical 



reentrant behavior, that is common to statistical systems [32l. I33l. I34ll35| . whenever com- 
peting interactions are present. Beyond the marginal value pjf ~ 7 .1A{3\E\'^ /SttG^ M^) 
there can always be found equilibrium states in dS case. 



In region I of Figure IHl there exist series 1 of equilibria of Figure 6(a) and in region 
II all equilibria of Figure |6(b)[ In the small upper gray shaded region there are the rest 



series of Figure 6(a) These type of equilibria exist only for values of the cosmological 
constant greater than a minimum value Pa""- This is the smallest value for which the 
cosmological force can keep marginally the whole matter at the edge. It can easily be 
calculated by equating the forces at the edge, assuming all matter is concentrated at R: 

= GM-npT'^R^ pX"" = P/4 (70) 

where p is the mean density. 

7.2. Comparison with Schwartzs child- dS space 

The two critical radii in Figure |9] resemble the two horizons of Schwartzschild-dS 
space, where the role of the cosmological horizon is played by Ria- The Schwartzschild- 
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Figure 9: The critical radius versus p/^ for fixed E, M in the microcanonical ensemble. There exist no 
equilibria in the unshaded region. With Rh is denoted the radius of the homogeneous solution. See text 
for details. 



dS metric can be written as: 

, 2 / 2GA/ SttG 2\ , 9 2GAf SttG 2\ "\ 2 2 ,^ n 

ds^--[l-^~ ^P.r^) dt^+[^~^- ^P.r^) dr^ + r^dn, (71) 

where pA is the 'mass' density of the cosmological constant (the energy density is k = 
p\c^ = ^q)- This metric has two horizons for p\ > 0; the black hole horizon Rbh and 
the cosmological horizon Re- Both are defined as the real roots Rh of the third order 
polynomial equation 

^ 2GM SttG . 

which, for various values of the cosmological constant p\, are plotted in Figure ITOl The 
resemblance with Figure [3] is too much striking to be considered accidental! It seems as 
though the reentrant phenomenon of Figure ^ is the closest Newtonian analogue of the 
horizons (Figure [TU)) in Schwartzschild-dS space. 

However, there is a big deference between the Schwartzschild-dS space and the New- 
tonian reentrant phenomenon. It is the opposite sense of the inequality for the instability, 
i.e. the stable region in the Newtonian case corresponds to the unstable region (i? < Rbh 
and R > Rc) of the Schwartzschild-dS space. 
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Figure 10: The horizons of the Schwartzschild-dS space versus the cosmological constant for a fixed 
mass. The horizon radius Rjj is measured in units of the Schwartzschild radius = For a given 

cosmological constant there are two horizons: the black hole horizon Rbh and the cosmological horizon 
Rc- The similarity with Figure |9] is striking! 



8. Canonical ensemble 




Figure 11: In dS, the dimensionless temperature versus the density contrast for fixed M, _R, pA in the 
canonical ensemble. Only the upper series (1 and 4) correspond to p(r) monotonically changing. At 
points B, except B4, an instability sets in. See text for details. 
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Figure 12: The critical density contrast versus pA fof fixed Af, K in the canonical ensemble, where p is 
the mean density of matter. An instability sets in for pi^/pn > {po/ pR)cr (except when {po/ pRjcr < 1), 
where a clumping phase transition occurs. 



In the canonical ensemble [l3| the instability sets in at different points than in the 
microcanonical ensemble. The series of equilibria expressed as function of temperature 
with respect to the density contrast can be seen in Figure [11] in case of dS. In AdS, just 
like flat case, there exist only one series of equilibria for a fixed cosmological constant, 



that is similar in form to series 1 of Figure 11(a) We note that Figure [TT] of temperature 
is very similar to Figure IH] of the energy. At points B, except i?4, of Figure [11] an 
instability sets in. Series AiBi and A4B4 are stable. A safe conclusion on the stability 
of the green series A2B2, A3B3, A5B5 and AqBq could not be reached, although the 
positive specific heat indicates stability, since we work in the canonical ensemble. The 
rest series are unstable. 



In Figure [T^ we can see that the critical density contrast for series 1 of Figure 11(a) 



decreasing for increasing cosmological constant, just like in the microcanonical ensemble. 
However, the instability sets in earlier in the canonical ensemble. It is well known in the 
flat case, and this fact is true with a cosmological constant, too, that the instability in 
the canonical ensemble sets in when the specific heat becomes negative. This negative 
specific heat region is stable in the microcanonical ensemble, while the instability sets in 
when the negative specific heat becomes positive again, in this ensemble (see Q])- 

The canonical ensemble is completely different then the microcanonical ensemble as 
far as the stability of the system is concerned. We alway consider the mass to be fixed 
in both ensembles. As can be seen in Figure I13[ the critical radius is decreasing with 
increasing cosmological constant and the region of the instability is now for radii smaller 
(and not bigger) then the critical radius R < Rcr- The critical radius is decreasing for 
increasing pA, because for an increase in pA the temperature is decreased, so that one 
should compress the system to balance out this destabilizing temperature decrease. The 
region of the instability changes, because in a compression although the pressure gradient 
is increased, it is weakened compared to the microcanonical ensemble, due to the heat 
transferred to the heat bath. In AdS case there is a marginal value of the cosmological 



24 




Figure 13: The critical radius versus pA for fixed /9, AI in the canonical ensemble. There exist no 
equilibria in the unshaded region. With Ru is denoted the radius of the homogeneous solution. See text 
for details. 



constant p^'^^ — —12 .32(3 / 8nG^ JVP (3^) beyond which no equihbrium states are possible. 

The critical temperature is decreasing with increasing cosmological constant, as can 
be seen in Figure [141 Increase of the cosmological constant acts as a stabilizer on the 
system due to the increase of the repelling force (or the decrease of the attracting force in 
case of AdS), enabling the system to be stable at lower temperatures. As the cosmological 
constant increases, it reaches a value for which the system can marginally be in static 
dynamical equilibrium. At this state all matter is still, i.e. T = 0, and is concentrated at 
the edge. This is point A in Figure [TH We have calculated this point earlier in equation 
(l70)) and found p™*" = p/4. For greater values the outward pointing cosmological force is 
increasing, enabling the mass to approach towards the center and to greater temperatures. 
The system undergoes a reentrant phase transition in the canonical ensemble. For a fixed 
cosmological constant at this region (after point A), there are metastable states for low 
temperatures up to some maximum critical value Ti, where a collapse phase transition 
occurs. For greater temperature, there exist no metastable states and the system suffers 
isothermal collapse. This happens up to some second critical temperature T2. For even 
greater temperatures, the equilibria are restored. 
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Figure 14: The critical temperature versus p\ for fixed M, R in the canonical ensemble, where p is the 
mean density of matter. In the unshaded region there exist no equilibria. This behavior indicates a 
reentrant collapse phase transition. 



9. Conclusions 

In the presence of the cosmological constant we find that, in the microcanonical en- 
semble, the critical radius, namely the Antonov radius, above which the system becomes 
unstable is increasing with increasing (negative or positive) cosmological constant, while 
in the canonical ensemble it is decreasing. The critical energy and temperature (less than 
which the system becomes unstable), as well as the critical density contrast {po/ PRjcr, 
are decreasing in the two ensembles. 

In dS case a new phenomenon is discovered, namely a reentrant phase transition; in 
microcanonical ensemble there emerges a second critical radius above which metastable 
states are restored, while in the canonical there appears a second critical low temperature, 
lower than which equilibria are restored. We stress out the similarity of the behavior, we 
have discovered, of the two critical radii in dS case in the microcanonical ensemble with 
the two horizons of relativistic Sschwartzscild-dS space. It seems that our 'dS case' is a 
Newtonian analogue of Schwartzschild-dS system. 

Another interesting feature of dS is the turning point of stability for the homogeneous 
solution, which resembles the Einstein's static universe, and the infinite numbers of non- 
uniform solutions at the homogeneous radius which suffer a transition from stability to 
instability when passing from solutions with density contrast po/pR < 1 to the ones 
with pq/ pn > 1. In addition, there exist multiple series of equilibria for a given positive 
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cosmological constant and fixed radius and mass. 

We stress out, that in this work the non-equivalence of ensembles in gravitating sys- 
tems is confirmed in the most dramatic way. For a fixed mass, in the microcanonical 
ensemble the instability occurs for radii larger than a critical value, while in the canon- 
ical the instability occurs for radii smaller than a critical value. That is because in the 
canonical ensemble, the pressure gradient that balances gravitation, is not drastically 
increased during a compression of the system, due to the heat transfer to the heat bath. 
In contrast, in the microcanonical ensemble where there is no energy loss, the pressure 
gradient is increased during a compression, drastically enough to balance gravity. 

Regarding applications to the physical world (positive cosmological constant), we 
state two interesting issues raised by our work. In the microcanonical ensemble, our 
analysis could have implications to the evolution of galaxy clusters, since many of them 
present a core in their centre and others a supermassive black hole. A quick calculation 
shows that, for the observed value of pA, the relevant to our stability analysis, dimension- 
less quantities 2pa/p and ii^G^ p\/?>\E\^ are of order unity for some typical values of 
regular galaxy clusters fl?]. This implies that the cosmological constant could influence 
the onset of the instability of galaxy clusters. In the canonical ensemble, the cosmo- 
logical constant could have an effect on the fractal structure of the Universe, since it is 
connected with the secondary instabilities 5] . Our analysis on the asymptotic behavior 
of the equilibria (section [S]), can be used to calculate the secondary instabilities. 

Let us close, noting that, even though there might be phenomenological connection 
with the physical world, our work is mainly focused and wishes to contribute on the 
theoretical understanding of the impact of an arbitrary cosmological constant term to 
the stability of gravitational systems. 

Appendix A. Antonov's proof that global entropy maxima do not exist 




Figure A. 15: The cylindrical coordinates for the second sphere. The first sphere, which is not drawn, 
has its center at the origin. 
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Let < a < 1 and aM be the mass uniformly distributed inside a sphere with radius 
ri and (1 — a)M the mass uniformly distributed inside a second sphere with radius r2. 
The distance of the spheres' centers is ri2. The distribution functions are constants 

^ _ aM _ aM f - ^ "-^^^ 

where vi, V2 the corresponding to the two spheres velocity bounds and the velocities v 

are equally probable by construction. 

Let us calculate the entropy of the system: 



S= - /caAflog/i - fc(l - a)log/2 = -fcM{aloga + (1 - a)log(l - a)} 

9M 
19^ 

The Newtonian dynamical energy Ui of the first sphere is 



+ 3A:M{(1 - a)logt;2r2 +aloguiri} - fcMlog:^^ (A.l) 



U,^- r G^^dm = -GM'a'^ C r'dr = -^GM^ — 
Jo r n Jo 5 ri 



and similarly for the second sphere 

f/2=-W(i^ 

5 r2 

Let our coordinate system be centered at the center of the first sphere and the center of 
the second sphere be at z = ri2. The cosmological dynamical energy for the first sphere 
is 

Uai = / pcpAU r = ^PA 4 — ^ / 47rr dr = --ttG p^Mar-^ 

Jo 3 |7rrJ Jo 5 

For the second sphere we use the cylindrical coordinates {T,(j),z) (see Figure lA.lSI) . Its 
cosmological dynamical energy is 

Jv 3 37rr| Jo Jri2-r2 Uo J 



-AnGpAMil - a) \ il(r?2 - (ri2 - r^f) + \rl 

The Newtonian dynamical energy between the two spheres is 

M2a(l_a) 

Ul2 — — Ct 

The kinetic energy of the first sphere is 

2 J 2 Jo Jo ^TT^rfvf 10 

and similarly for the second sphere 

i^2 = |:(l-a)Mi;2 
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Summing it all up, the total energy is 



E 



-GAP i — + 



5 

AttG 



n 



MpA 



r2 
(1 



5 a(l — a) 
3 



ri2 



a) 



5 1 

6 r2 



+ {avf + {1 - a)vl} 

We keep constant the quantities ri , V2 and L 



1 r 10 

a\3M 



E 



GA''Pa(l - a) 



ri2 



(^?2 - {ri2 - ra)') + -r^ 

(A.2) 

(1 — a) log r2. Equation (jA.2[) gives 
(1 - a)vi + 



2GM— 



SttG 



PA ar^ + {l-a) 



5 1.3 . 



r2 



2GL^M 

''2 (log U2f2)^ 

(A.3) 



For r2 — > equation (|A.2p gives a — > 1 in order to keep the energy finite. For r2 — > 0, 
a — >■ 1 equation (jA.3p gives ui oo. For r2 — > 0, a — > 1 and -ui — > oo equation (jA.l[) 
gives S" — >■ oo. Therefore there exists a configuration with finite energy for which entropy 
is not bounded from above, or equivalently, matter can always be redistributed in such 
a way keeping energy fixed and increasing the entropy. 



Appendix B. Poisson equation with A 

Let us calculate the equation of the gravitational potential (j) in the presence of a 
cosmological constant A, in the Newtonian limit. The Einstein's equations are 

R^, ~ ^RSi: - A<5^; = ^T!; (B.l) 
In the non-relativistic limit, the energy- momentum tensor is 

T^' ~ pc'5^51 

Contracting the Einstein's equations we get R = —SttGT — 4A and substituting again in 

= l^j/. _ (B.2) 



equation (jB.ll) we get 



In the weak field limit only the time components of the Einstein's equation survive. For 
slowly moving particles it is 

^2 i 1 



dt^ - - -00^-00- ^2- 

and in the static weak field limit it is 

^0 = -^00 — -j^dad" goo = ^ctFqq = ^V^<i 



Then, the time-time component of equation (jB.2l) gives 

V'^cj) = 47rGp - SttGpa 

where 

Ac2 

PA 



SttG 
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Appendix C. Derivation of ST 

Let us give a useful expression for the dynamical energy 

^ = \j P4'Nd^r + j p(j)Kd^f^ i j p{(f)- (j)A)d^f+ J p(t>Ad^f 

= \ j pc^d'r + i I pc^Ad^'r (C.l) 
We assume pA to be fixed. Using (jC.l|) we have 



5E = 5K + 5U=—5T+{5p)— + {54>)— + -[{5p)— + {54>)—] C/ + 0(3) 
= -rrST + - / d^r{pd<pN + (t>NSp + (f>A5p + + t; I d^" 



= + 1 dV- (^0<5p + ]^5p5<P^ (C.2) 

In the third raw we used the identity 



d^r{(t>N5p + p5(t)N) = / dV(5p(r) / d^r' 4, + / dVp(r) / 



"^Jd^^J d^f'6p{r)-^^^ = J Sr2(j)N5p (C.3) 



The constraint 5E — gives by use of (IC.2P : 

Appendix D. How to determine the unstable branch near a turning point 

We search for a solution of the problem ((36)) for ^ = 0. Let T be the operator 

f d ( 1 d\ , G 



dr \47rpr^ dr J Tr"^ 

Let Vt be the trial value we use to solve If it does not correspond to the solution 

-Fo, it would correspond to some other solution Fn for a different eigenvalue 

TFn ^ ,^Vt + inFn (D.l) 
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so that 



Vt= f dr<j>'F^ 
Jo 

Let V be the value of the integral 



V= I dr(t)'Fo 



where Fq is the solution of the problem 



TFo = -^^Vt (D.2) 



_Fo will correspond to a solution with a zero eigenvalue only if 

V = Vt 



Equation (jD.2p gives 

r-fl OT/_ rR 



and 

3A/r2 7o ^ " 3Mr2 



dr F^TFo = ttI^ / dr = 77tI;;^^t 



Let us prove that 



drFnTFo = / drFoTF^ (D.3) 



We have 



R rR 



drF^TFo - / drF,A^{-^^]+^Fo 







c?r \47rpr2 dr / Tr^ 



1 p dFo 



Anpr^ dr 
1 di^. 



[ dr \4TTpr^ dr J Tr'^ " ^ 



47rpr2 dr 



"'0 
R nR 

■Fn + / drFoTF^^ / drFoTF^ 

"'0 



Using equation (jD.3l) and subtracting (jP.ip from (ID.2P we get 

2V^ ( V \ 

TFo - TF^ = ^ ^ (^^ - ij - -6. ^^oF„ 

Near the turning point it is /J* dr F^Fn — dr Fq so that an instability > 0) sets 
in when V < Vr- 
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Appendix E. Calculation of the total energy 



Let us calculate the expression for the total energy with no use of the virial theorem. 
First, we need the value (/)(0) at the origin. Equation (jlOp gives 

0(0) = -G / dr I^Anr^ = -G / ^ poe^^47r^=== ^ = - / dxxe'y 







r Jo ^/AnGpoP ^/AnGpo^ 

In the followings we use equation (jC.l[) . We have 



E 




ER 


3 




2B 




3 




W 




3 




2B 



2/3 2]'^^ 2 
GAP 2 y ^^^^ + ^^^0))'^- GM^2— io ^'^^^^ 

i^)2R^4^FG^2j, ' +2^f^ 12(GAf/3)2 po io ' 

1 /"^ A 1 /"^ 

dxx^e^^iy x'^) / dx xe~^ (E 

6 ' 2B.L ^ 
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